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The process in which a photon splits into three after the collision with a free electron (triple 
Compton effect) is the most basic process for the generation of a high-energy multi-particle entangled 
state composed out of elementary quanta. The cross section of the process is evaluated in two 
experimentally realizable situations, one employing gamma photons and stationary electrons, and 
the other using keV photons and GeV electrons of an x-ray free electron laser. For the first case, 
our calculation is in agreement with the only available measurement of the differential cross section 
for the process under study. Our estimates indicate that the process should be readily measurable 
also in the second case. We quantify the polarization entanglement in the final state by a recently 
proposed multi-particle entanglement measure. 
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Introduction. — The triple Compton effect is a funda- 
mental process of light-matter interaction, in which a 
photon splits into three upon collision with an electron, 

c~+7 — > c~ +7 + 7 + 7. (1) 

The cross section of this quantum electrodynamical 
(QED) process is of fourth order in the fine-structure 
constant a w 1/137.036. Currently, there is no general 
theoretical treatment of the triple Compton effect in the 
literature. In principle, there is no limit of the number of 
photons that can be coherently emitted when an electron 
interacts with a photon, but the only processes which 
have been discussed so far are the usual, second-order 
(in a) single Compton effect [1] and the third-order dou- 
ble Compton effect [2]. Only the single Compton effect 
has a classical limit. To the contrary, both the double 
and the triple Compton effects are intrinsically quantum 
processes which cannot be described by classical electro- 
dynamics. With strong lasers becoming available, the 
nonlinear generalization of Compton scattering in which 
several laser photons are absorbed have been vigorously 
discussed; see [3, 4] for recent investigations of the non- 
linear Compton effect, and [5-7] for a discussion on the 
nonlinear double Compton effect, where two photons are 
coherently emitted. 

The three photons emitted in the process (1) originate 
from the same initial photon, are emitted at the same 
time, and are therefore quantum mechanically entangled. 
The creation of an entangled state of three photons is an 
important goal in quantum information. The conven- 
tional way of generating entangled triple-photon states 
is by employing nonlinear crystals [8-10], but one can 
imagine other sources, not yet experimentally realized, 
such as electron-positron annihilation [11, 12]. In this 
Letter, we propose the triple Compton effect as an al- 
ternative source of entangled photon triplets, which by 
suitable optimization could be competitive both in terms 
of production rate and the degree of entanglement. Espe- 



cially interesting is the possibility of creating correlated 
photons with high energy in the GeV range, via triple 
Compton backscattering on a relativistic electron beam. 
We will show that a advantageous setup for such an ex- 
periment is an x-ray free electron laser (XFEL), provid- 
ing a high-flux, high-energy photon beam together with 
a GeV electron beam. 

The only previous experimental study of the triple 
Compton effect that we are aware of is Ref. [13], where 
the differential cross section was measured for a well- 
defined interval of the solid angle, as defined by the de- 
tectors which were arranged in a symmetric angular con- 
figuration. On the theoretical side, the only preceding 
investigation is [14], where the total cross section in the 
limit of ultra-relativistic initial photon energy was stud- 
ied. In our treatment of the problem, we can evaluate 
the differential cross section for arbitrary values of the 
directions, energies and polarizations of the emitted pho- 
tons. On the other hand, the double Compton effect is 
rather well studied, both experimentally [15-18] and the- 
oretically [2, 19, 20]. In addition, several other processes 
have been studied where two photons are produced in 
the final state, such as double bremsstrahlung [21-23], 
bound state decay [24, 25], and laser- induced photon 
splitting [26], but comparatively little is known about 
triple photon production. 

Theoretical formulation. — Unless stated otherwise, 
natural units with h = c = 1 are used, and m denotes the 
electron mass. Four- vector products are denoted with a 
dot (i.e., we have a ■ b = a^b^ = a b° — a ■ b for two 
four- vectors a and b). The contraction with the Dirac 
matrices 7^ is written with a hat, a = a°7° — 0-7. 

We label the four-momenta of the incoming electron 
and photon with pi = {Ei,pi) and fco = wo^o = 
wo (1, no), respectively, and the four- momenta of the out- 
going electron and photons with pj = (Ef,pj), and 
kj = ujjnj = u>j(l,rij), j = 1,2,3, respectively. When 
explicitly evaluating the cross section, we take the vectors 
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FIG. 1. One of the 4! = 24 contributing Feynman diagrams of 
the triple Compton effect. To obtain the total amplitude, we 
have to add the diagrams corresponding to the permutations 
of the insertion order of the photons k , fci, k 2 , and fc 3 . 



kj = ujjTij in the lab frame, in spherical coordinates with 
the polar axis directed along fco, i.e., fco = u>o(0, 0, 1) and 
kj = ujj (sin 9j cos <j>j , sin 9j sin <j)j , cos 9j ) . The amplitude 
M for the triple Compton effect is formed from a coher- 
ent sum of 4! = 24 Feynman diagrams, one of which is 
shown in Fig. 1. 
We have 

M = m 2^4, (P/h e ?(3) 727^3^2 e «(2) ^27 



g?(0 + m 



X ^ (1) g 2 (0 - m 2 g «0)"ri(Pi) . 



(2) 



where the sum runs over all the 24 permutations £ 
of (0,1,2,3), and r/ is the final electron spin. The 
momenta q n entering the propagators are defined as 

« n (0 = Pi + E"=o(-l) 5o£0)+1 %a), where is Kro- 
nccker's delta function. (One adds to pi the four- 
momentum of the incoming photon fco and subtracts 
the four-momenta of the outgoing photons, according to 
the relevant permutation £.) In the expression for the 
amplitude (2), we have introduced the positive-energy 
bispinor u r (p), with the spin index r = 1,2, which we 
use in the conventions of Chap. 2 of [27] (the normal- 
ization is u r (p)u r (p) = u I (p) 7° u r (p) = 1). The po- 
larization four- vectors ej = (0, Cj) satisfy ej ■ kj =0 
for j = 0,1,2,3. As an explicit set of basis vectors 
for ej , we take e] = (cos 9j cos <fij , cos 9j sin <f>j , — sin 9j ) 
and e 2 = (— sin <f>j, cos <pj, 0). The differential cross sec- 
tion follows from the usual rules of QED [28] as 



da 

dwidw 2 df2id£l 2 d0 3 
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a 4 1 cjiw 2 w 3 |M| 2 
(2tt)4 rr^EfPi ■ k \K\ 
xe(cj 3 )B{E f -m), (3) 



where O(-) is the step function, dfij = sin9jd9jd<pj is 
the differential solid angle of photon j, and we have in- 
troduced (T5 as an abbreviation of the five-fold differential 
cross section. In Eq. (3), the final four-momentum of the 
electron is fixed by energy-momentum conservation as 
Pf = pi + fco — fci — &2 — k^, and the energy of photon 



w 3 



Pi ■ (fc - fci - fc 2 ) - fcp ■ (fci + fc 2 ) + fci • fc 2 
n 3 ■ {Pi + fc - fci - k 2 ) 



(4) 



The factor K in (3), arising from the final delta function 
integration over lo 3 reads K = 1 + ^-(fci+fca-fco-pO+h* _ 

The differential cross section (3) depends on eight con- 
tinuous variables, two angles 9j, <fij for each emitted pho- 
ton and the energies uj\ , lo 2 of two of the emitted photons. 
The energy L03 of the third photon is restricted by the 
kinematical constraints and is calculated by Eq. (4). In 
addition, we have six more discrete variables which can 
take any of two values, namely, the polarization vectors 
6j of the photons and the spins r^j of the electrons. 

A couple of general remarks about the expression (3) 
are the following. Similarly to the double Compton ef- 
fect [28], the cross section vanishes when all three pho- 
tons are emitted parallel to the incoming photon. In the 
rest frame of the electron, this implies u>i + cj 2 + UI3 = ujq 
according to Eq. (4), and therefore corresponds to the in- 
coming photon splitting into three photons without any 
interaction, which is not allowed. Whenever either of 
u>i, 0J2 or CJ3 goes to zero, 0-5 diverges as l/wi, 2 .3, while 
the radiated energy remains finite in the infrared. This 
is the well-known infrared catastrophe of QED, which 
can be cured by adding radiative corrections [29]. In 
the current case, the divergences at uj n — >• would can- 
cel against radiative corrections to the single and double 
Compton effect. The evaluation of such corrections is 
beyond the scope of this Letter. We will calculate the 
differential cross section sufficiently far from the infrared 
divergences, corresponding to a specific experimental de- 
tector threshold. Our calculations therefore have a rela- 
tive accuracy of the order of the fine-structure constant. 

Evaluation of the differential cross section. — The eval- 
uation of 0-5 is performed numerically, by employing an 
explicit representation of 7^ and u r (p). This approach 
is advantageous [30] if one is interested in polarization- 
resolved cross sections, in which case the analytic eval- 
uation of cr 5 does not simplify. We calculate the differ- 
ential cross section (3) for two different experimental se- 
tups. The first is the same as considered in Ref. [13], 
were a measurement of the triple Compton cross sec- 
tion is described. In this setup, a gamma photon of en- 
ergy ui = 0.662 MeV impinges on an electron at rest, 
Ei = m, and photons are detected in coincidence at 
9i = 9 2 = 9 3 = it/2, and (f>i = 2tt/3, (f> 2 = 4tt/3, <fe = 0. 
The photon energy threshold was u n > e — 13 keV. 
What was actually measured was the differential cross 
section averaged over the solid angles subtended by the 
detectors, 

^ = ^\ / dnidJ^dfi^ /"dwidw2<75, (5) 

spin, pol.-'^ 1 J 

where the solid angle is 17 = 0.378 sr, and the energy 
integration is over wi j2 .3 > e. Performing the integration 



in (5) over the angular intervals A(/)j = <j>j±y/Q/2, A6j — 
6* J ±arcsin(v / ^/2) so that j Ae dcos^- f A ^ d<//- = O, and 
over photon energies greater than e, we obtain 

(a) =4.1 x 10~ 9 b/sr 3 . (6) 

This value should be compared to the measured value of 

(a) oxp = (8.1 ± 2.4) x 1(T 9 b/sr 3 , (7) 

which includes an experimental uncertainty of 30% [13]. 
Although the calculated value lies outside the estimated 
experimental error bar (by 1.6 standard deviations), we 
believe that our calculation rather can be regarded as 
a confirmation of the measurement in Ref. [13], given 
the utmost difficult nature of the experiment. We have 
extended the analysis relevant to the experiment [13] to 
the energy interval w € (10~ 2 , 10 2 ) MeV; details will be 
presented elsewhere. 

The second example we consider is a triple Compton 
backscattering scheme: The head-on collision of a rela- 
tivistic electron, at an incident energy of Ei — 5 GeV 
and Pi = (0, 0,—y/ Ef — to 2 ), with an incoming photon 
of energy uiq = 1 kcV. Such parameters are realized in 
an XFEL, for example at the LCLS in Stanford [31], and 
would require the reflection of the x-ray beam to collide 
with the electron beam from the accelerator. In this sit- 
uation, the photons are backscattered and emitted in a 
narrow cone around the propagation direction of the elec- 
tron, 6j ~ 7T — Ei/m. We estimate the total cross section 
as cr^c — 2 x 10~ 5 b, compared to a^c = 1 x 10~ 3 b and 
CTgp = 3x 10 -2 b, where for DC and TC, we have assumed 
that all photons with energies larger than the assumed 
photon energy threshold e = 50 MeV are detected. If we 
adopt pulse parameters available at the LCLS [31], i.e., 
2 x 10 13 photons per pulse, 10 9 electrons per bunch, a 
transverse bunch size of 40 fim, perfect transverse over- 
lap of the two pulses, and a repetition rate of 120 Hz, we 
obtain 3 triple photon events per second. 

One example of the fully differential cross section in 
the XFEL setup is shown in Fig. 2. In this figure, we 
assume a configuration centered near the incoming elec- 
tron axis, with 9j = it — 1.5 x 10~ 3 and 4>j = 2j7r/3, 
j = 1,2,3. We evaluate 05 as a function of io\ and o>2- 
We sum over the electron spins, but assume a linearly 
polarized XFEL beam, and also fix the polarization of 
the emitted photons. We show 05 only for 0^2,3 > s. 
For photon energies smaller than e, the differential cross 
section is set to zero. It becomes clear from Fig. 2 that 
high energy photons in the GeV range can be achieved 
in this setup, which is interesting since it could provide 
correlated photon triplet states, with quantifiable entan- 
glement as shown below, in an energy range which is 
far beyond what can be produced with down-conversion 
of photons in a nonlinear crystal [8-10]. The rate of 3 
events/s obtained above is not small compared to the 
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FIG. 2. (color online). The five- fold differential cross sec- 
tion at 9j = it — 1.5 x 10~ 3 and <j>j = 2jiv/3, j = 1, 2, 3 as a 
function of o?i and CJ2 . The color indicates the value of the 
cross section log 10 (| ]TJ spin os/bMeV -2 sr -3 ) summed over 
the electron spins. We have Ei — 5 GeV and Wo = 1 keV, 
and the incoming photon beam is linearly polarized in the re- 
direction. The different panels correspond to different polar- 
ization of the final photons, we have (ei, £2, £3) = (e^, e*, 63), 
with jkl = (a) 111, (b) 211, (c) 121, (d) 112, (e) 221, (f) 212, 
(g) 122 and (h) 222. The black dashed line in panel (a) shows 
the curve along which the energy of photon three is at the 
assumed detector threshold, 0J3 — e — 50 MeV. 

experimental results in [9], where photon triplets pro- 
duced from nonlinear down-conversion were detected at 
an event rate of 5 per hour. The polarization of GeV pho- 
tons can be measured using coherent electron-positron 
pair production in aligned crystals [32]. 

Entangled photon triplets. — The three photons in the 
final state are simultaneously emitted and their polariza- 
tions are inevitably entangled due to the quantum nature 
of the process. In general, the mixed polarization state 
resulting after the interaction can be described by the 
density matrix g, which has matrix elements 

(AiA 2 A 3 |g|A' 1 A 2 A^) =NJ2 M(AiA 2 A 3 ) M*(\[\' 2 \' 3 ), 

n,r f 

(8) 

where we have written M(AiA2A 3 ) = M(ei = £^,62 — 
e2 2 ,€3 = Eg 3 ), Xj <E {1,2}, and the normalization con- 
stant N is chosen so that Tr(p) = 1 (we sum over initial 
and final electron spins and rt). It is still an open 
problem and subject of active research how to quantify 
the entanglement of a given multi-particle state g (see 
Refs. [33-35]). In the current study, we employ the entan- 
glement measure put forward in [35] in order to estimate 
the degree of polarization entanglement present among 
the final three photons. Briefly, an entanglement witness 
W is found by minimizing the trace Tr(Wg) with W such 
that for all subsets s,W — P s +Qj s , where T s denotes the 
partial transpose with respect to the subset s [36], and 
the matrices P s , Q s , 1 — P s and 1 — Q s should have pos- 
itive eigenvalues. Then, t(q) = — Ti(Wg) is a measure 
of the tripartite polarization entanglement present in g. 
In particular, r = for states which are not genuinely 
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regime is possible based on the triple-Compton effect as a 
basic three-quanta emission process. The triplet photon 
final state is polarization entangled. 

We acknowledge support from the National Science 
Foundation and from the National Institute of Standards 
and Technology (NIST precision measurement grant). 



FIG. 3. (color online). Entanglement measure r as a function 
of the photon energies oji and wj, for (a) ujq = 0.662 MeV, 
E % = m, 6j = tt/2 ~ F <t>i = 2j7r/3, j = 1,2,3, e = 13 keV, 
and (b) w = 1 keV, E t = 5 GeV, 6j = tt - 1.5 X 10~ 3 , 
0j = 2j7r/3, j = 1,2,3, e = 50 MeV. In both cases spin was 
summed over, and the initial photon beam was assumed to 
be linearly polarized in the x direction. The differential cross 
section, and thus also r, was set to zero for photon energies 
below the detector threshold e. 



tripartite entangled. As an example, t(^ghz) = 1/2 for 
the Greenberger-Horne-Zeilinger (GHZ) state £>ghz — 
|GHZ)(GHZ|, with |GHZ) = (|111) + |222))/-\/2 [37], so 
that r = 1/2 may be regarded as maximum entangle- 
ment. 

In Fig. 3, we show the value of r, evaluated as a func- 
tion of uii and uj^ at (a) ljq — 0.662 MeV, Ei = to, 
#1,2,3 = 7r /2 — 1, and <j)j = 2j7r/3, and (b) w = 1 keV, 
e[ = 5 GeV, 6^2,3 = tt-1.5 x 1CT 3 , and ^ = 2j7r/3 (the 
same situation as in Fig. 2). The numerical evaluation of 
r was performed with PPTMIXER [35], available at [38]. 
In both cases a non-zero value of r for almost all values 
of uj\, uj2 where 05 > shows that the three photons 
are indeed entangled. Almost maximum entanglement 
t = 1/2 is achieved at certain places in the uj\U)2 plane. 

Conclusions. — We have presented a theoretical study 
of the triple Compton effect, where a photon is split into 
three after colliding with an electron. Our formulation of 
the problem enables us to evaluate the differential cross 
section at arbitrary angles and energies of both initial and 
final particles. We verify the 44-year old experimental 
result reported in Ref. [13]. Theory and experiment are 
not in perfect agreement, but the discrepancy of roughly 
1.6 standard deviations is not significant [see Eqs. (6) 
and (7)]. Additional measurements are needed to clar- 
ify this issue. A straightforward generalization of the 
formalism to a Compton backscattering geometry then 
leads to theoretical predictions for the triple Compton 
effect for typical parameters at modern XFEL facilities, 
as shown in Fig. 2. Finally, while it is intuitively rather 
clear what a two-particle entangled quantum state is, the 
quantification of three-particle entanglement is a much 
more subtle problem. The entanglement measure r has 
been proposed in Ref. [35], and corresponding theoret- 
ical predictions in Fig. 3 show that the generation of a 
multi-particle entangled photon state in the MeV or GeV 
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